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Abstract 

We explore the possibility of creating of solitons in gravitating gas. 
It is shown that the virial arguments does not put an obstacle for the 
existence of localized static solutions. The simplest toroidal soliton of 
gravitating gas could be the explanation of the peculiar galaxy named 
Hoag’s object. 


1 Introduction 


In the recent paper [1] we have considered the Hamiltonian formalism for 
fluid and gas based on the Lagrangian description. It was pointed out in this 
paper that apart from other advantages, the Lagrangian description, which 
uses the trajectories of the particles of fluid (gas) as the dynamical variables, 
is the most convenient for the introduction of interaction. In particular it was 
demonstrated in [1] how the introduction of the electromagnetic interaction 
of particles which constitute the fluid, provide us with the theory of plasma. 
In the present paper we are going to consider in analogous way the theory of 
gas of particles which interact with each other through gravitational Newton 
potential. The system of such particles could be considered as a model for 
the motion of stars in a galaxy when the gravitation interaction prevails all 
other interaction. The total number of stars in typical galaxy is of the order 
of 10^^ —10^^, so it may be reasonable to consider this collection of’’particles” 
as a gas. 

The simplest model, which is usually used for numerical simulation of 
N-body model of galaxy is described by the Hamiltonian 


H = 


^ V- 

E- 

i 


N 

-7E' 




where Xi,Pi are canonical coordinates of particles (stars) with masses m*. 
The model we are going to consider is based on the dynamics described by 
H with the assumption = m, when N —>■ oo. For this limit there appears 
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a natural desire to consider a continuos distribution of the particles as it is 
done in the theory of fluid or gas [2], 

Apparently, the attractive potential provided by gravity will produce es¬ 
sential difference of the properties of such a media. The potential part of 
the energy will produce a collapse, which will be prevented by kinetic term. 
Such a competition sometimes results in a creation of the steady solution — 
a soliton. 

It is important to mention here that the solitons in gas and fluid is in¬ 
trinsically different with e.g. single solitons considered in particle physics. 
In the case of fluid or gas only the held of density and velocities are time 
independent, while the constituents of the media are moving. The simplest 
example of such a soliton is a tornado. 


2 Description of the Model 


In huid (gas) mechanics there are two different pictures of description. The 
hrst, usually refereed as Eulerian, uses as the coordinates the space dependent 
helds of velocity and density. The second, Lagrangian description, uses the 
coordinates of the particles x{^i,t) labeled by the set of the parameters 
which could be considered as the initial positions ^ = 0) and time 

t. These initial positions ^ as well, as the coordinates x{^i,t) belong to some 
domain D R^. In sequel we shall consider only conservative systems, 
where the paths of different particles do not cross, therefore it is clear that 
the functions x{^i,t) dehne a diffeomorphism of D C and the inverse 
functions ^{xi,t) should also exist. 


XjiCiR) 

^j{Xi,t) 


^ ^ = Xj, 

= ir 


The density of the particles in space at time t is 


p[xR) = J d^Cpo{Ci)S{x - x{Ci,t)), 


( 1 ) 


( 2 ) 


where Po(0 is the initial density at time f = 0. The velocity field n as a 
function of coordinates x and t is: 


v{Xi,t) = x{^{Xi,t),t), 


( 3 ) 
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where ^{x,t) is the inverse function (1). The velocity also could be written 
in the following form: 


v{xi, t) = 


f d^^po(^i)x(^i, t)6(x - x((i, t)) 

I d^^Po{Ci)S{x-x{Ci,t)) 


(4) 


or 


p{xi,t)v{xi,t) = J d^^po{Ci)x{Ci,t)S{x - x{^i,t)). (5) 

Let us calculate the time derivative of the density using its dehnition (2) : 

d 


p{xi,t) = J d^Cpo{^i)-^S{x - x{Ci,t)) 


^ J d^CpoiCi) (-^(6, f)) ^)) 

= I d^^Po{Ci)x{Ci,t)d{x - x{^i,t)) 
d 

= -l^P[Xi,t)v[Xi,t) 

ox 


( 6 ) 


In such a way we verify the continuity equation of fluid dynamics: 

p{xi, t) + d{p{xi, t)v{xi, t)) = 0. (7) 


Using the coordinates x{^i,t) as a conhgurational variables we can con¬ 
sider the simplest motion of the fluid described by the Lagrangian 

i = (8) 

The equations of motion which follow from (8) apparently are 

mx{^i,t)=0 (9) 

Now let us hud what does this equation mean for the density and velocity 
of the fluid. For that we shall differentiate both sides of (5) with respect to 
time 

d [ " 

—p{xi,t)v{xi,t) = J d^Cpo{^i)x{^i,t)S{x - x{^i,t)) 

+1 d^^po{^i)x{^i,t)^6{x - x{^i,t)) (10) 
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The first term in the r.h.s. of (10) vanishes due to the equations of motion(9) 
and transforming the second in the same way, as we did in ( 6 ) we arrive at 


d d 

t) + — t)v{xi, t)vk{xi, t)) = 0 

Let us rewrite (11) in the following form: 

d 


( 11 ) 


v{Xi, t) p{Xi, t) + — (^p{Xi, t)Vk{Xi, f)) 


d 


+p{xi, t) \y{xi, t) + Vk{xi, t)—v{xi, t) 


= 0 . 


( 12 ) 


The hrst term in (12) vanishes due to the continuity equation, while the 
second gives Euler’s equation in the case of the free flow: 


d 

v{xi, t) + Vk{xi, t)- — v{xi, t) = 0 
UXh 


( 13 ) 


In order to get the usual Euler equations for fluid or gas with the internal 
pressure we need to add to the Lagrangian ( 8 ) the ’’potential” part, as it 
has been shown in [1], This ’’potential” part describes the repulsive inter¬ 
action between particles which constitute the media (gas or fluid). As was 
announced in Introduction, we are going to introduce another interaction 
— gravitational attraction between particles. Eor that we have to modify 
the Lagrangian ( 8 ) in the following way: 


L 




mx 


7 f Po(^»)po(40 

2 7 ^ Nf(e*,f)-i^(e),t)r 


(14) 


where 7 denotes the gravitational constant. The equations of motion, which 
follow from the Lagrangian (14) have the form: 


mx{^i, t)+'y 


d^ePoiQ 


x{^i,t) -xXi,t) 
x{Cut) -xXi,t)\^ 


0 


(15) 


Translating equation (15) on to the language of Euler variables, as has been 
done above we arrive at the following set, including the continuity equation: 


v{Xi, t) + Vk{Xi, t)—v{Xi, t) 
Kxi,t) + ■^(^p{Xi,t)v{Xi,t)) 


mdx J \x — y\'’ 

0 . 


( 16 ) 

( 17 ) 
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Note that the r.h.s of the equation (16) in the case of ordinary gas or fluid 
is expressed through the internal pressure p{xi)] 


v{Xi,t) +Vk{Xi,t)- — v{xi,t) = - ———p{xi 

dxk p{xi) dx 


(18) 


The set (16),(17) dehnes the evolution of initial distribution of p{xi, to),v{xi, to) 
of gravitating gas and, besides it could be used to hnd the static conhguration 
of this gas for different bonndary conditions. In particular we can explore the 
possibility of the existence of the static isolated conhgurations of gravitating 
gas. Isolation here means that density p{xi) vanishes at inhnity. Note, that 
for usual fluid, gas or plasma these kind of solntions are forbidden dne to 
virial argnments, known in the case of plasma as Shafranov’s theorem [3],[4]. 
For the gas, describing by equation (18) this theorem could be proven as 
follows. First, using the continuity equation let us rewrite equation (18) for 
static case in the following form; 


A. 

dxk 


p{xi)vk{xi)vj{xj) + 5jkP{xi 


= 0 . 


Integrating (19) with Xj over we obtain: 

d 


0 = f (Rxx 


f (Rx 


A. 

dxk 


^ dxk 


x, 


p{xi)vk{xi)vj{xj) + djkp{xi 
p{xi)vk{xi)vj{xj) + 5jkP{xi 


J (Rx6jk p{xi)vk{xi)vj{xj) + SjkP{xi 


(19) 


( 20 ) 


For usual gases p{xi) ~ p'^{xi),'y > 0, therefore the integral over divergence 
will vanish for isolated solntions for which p{xi) 0 when |a;| —cx) and we 
the obtain the following equation: 


(fx 


p{xi)iP{xj) + 3p{xi) 


= 0 . 


( 21 ) 


Apparently, this eqnation could be satished only for the case p{xi) = 0, i.e. 
there is no isolated in the above formulated sense, static solntions of the 
equation (18). In order to have a static solntion of (18) we need to change 
the bonndary condition p{xi) —0 to the condition p{xi) —> pas , where pas 
is asymptotic nniform density. 
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Now we shall show that the arguments of this theorem bring no obstacles 
for gravitating gas. For this we again will rewrite the equations (16) for the 
static case, using continuity equation in the following form: 


- d 

.dxj 


p{Xi)Vk{Xi)Vj{Xi] 



0 , 


Integrating (22) with x^ over we obtain: 


( 22 ) 


/ 


d^xxk 


- d 

_dxj 


p{Xi)Vk{Xi)Vj{Xi 



0. (23) 


Consider the first term of the integrand in (23). Integrating by parts as above 
and taking into account the asymptotic conditions for p{xi) we obtain: 


/ 


d^xxf,-^(^p{xi)vk{xi)vj{xi)^ ^ - jd^xp{xi)if{xi 


( 24 ) 


Integration of the second term of the integrand in (23) is straightforward, 
yielding 


m 



7 

2m 


\x - y\ 


In such a way from equation (23) we obtain: 


(fxp{xi)iP{xi) — / cfxcfy 

^ 11L J 


p{Xi)p{Xi) 

\x-y\ 


= 0 . 


(25) 


(26) 


This relation apparently could be satisfied for non-trivial conhgurations of 
p{xi),v{xi). The energy functional, corresponding to the Lagrangian (14) 
has the following form: 


E = 


m 


d^xp{xi, t)iP{xi, t) 



|i-! 7 l 


( 27 ) 


The two terms of the energy functional have clear interpretation as kinetic T 
and potential U parts of energy and equation (23) expresses famous ’’virial 
theorem” [5]: 

2T = -U. (28) 
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Note here, that in the ’’virial theorem” equation (28) holds true for mean 
values of kinetic and potential energies, while in our case of static solutions 
there is no need to average over time. 

Using the relation (28) we can easily hnd the total energy of static con- 
hguration of the gravitating gas: 

Ecstatic 


So, the total energy of the static solution is negative, as for the bound state 
of Kepler problem. 

3 Properties of Static Solutions 

The equations which dehne our static conhgurations of gravitating gas have 
the following form: 


m 


d^xp{xi)v^{xi) = 


( 29 ) 


'I /■ J 3 .j 3 . 


(Txcry- 


\x-y\ 


Vk[Xi, 


A. 

dxi 


-ViXi) = 


7 5 /■ ,3„, PiVi) 


m dx 




x-y\’ 


A 

dx 


(^p{xi)v{xi)'j = 0 


Taking divergence of the first equation : 


d 


Xi 


A. 

dxi 


7 a /■ PiVi 


v(xi)) = —A / d-^y 
' m J 

7 


\x-y\ 


= —{-ATl)p{Xi), 

m 


( 30 ) 


( 31 ) 


we can write the whole set of the equations for static conhguration in a pure 
local form: 


d(vk{xi 


A_ 

dxi 


-V[Xi) = 


d 

— {^p{xi)v{xA) = 0, 


dx 


d 

dx (vk{xi)—v{xi)'^ = 0 , 


dx 


A 

-An—p{xi), 

m 


(32) 
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where the last equation requires the expression Vk{xi)-^v{xi) to be a gradi¬ 
ent. 

Now we are going to derive an important inequality, which bounds the 
{—Estatic) of any solution of (32) from below. For this let us introduce the 
notation for the potential U{xi): 


U{x,) = - [ (33) 

J \x-y\ 

Integrating by parts we obtain the following relation: 

[ d^x(dU{xi)f = An [ ( 34 ) 

J J uZ' y 

Using (34) we can write the expression for the energy given by (29) of any 
static conhguration, which is the solution of (32) in the following form: 


e.,.h, = J d^x[au(xt)f 


(35) 


Substituting into (35) the expression for the gradient of the potential U{xi) 
from the first equation (30), we obtain the expression for the Estatic only 
through the held of velocity: 


-E 


static 



(36) 


This form is most convenient for the derivation of desired inequality. Now let 
us consider a function f{xi) from the Hilbert space IU 2 ^, which consists of all 
measurable functions on R^, which have at least one derivative and square 
integrable on together with its derivatives. In particular we assume that 
the density p{xi) belongs to W 2 . Taking into account the equation (31) we 
have 


\Jd^xp{xi)f{xi)\ = J d^xf{xi)dk(vj{xi)djVk{xi)'^\ 

= I d^xdkf{xi)vj{xi)djVk{xi)\ 

< -^^{1 d^x{dkf{xi))^y^^(^J d^x{vj{xi)djVk{xi)fy^‘^, (37) 


where on the last step we have used Cauchy inequality. From (37) we imme¬ 
diately obtain the inequality: 


-E 


static 


(/ d^xp{xi)f{xi)y 

“ J d^x{dkf{xi)y 


(38) 



which is valid for any f{xi) from W 2 and is saturated for /(x*) = U{xi). 
Indeed, let us calculate the derivative of the functional in the r.h.s of (38) 
with respect to f{xi): 

6 (f d^xp(xi)f(xi)') 

df{xi) J d^x{dkf{xi)y 

^ {jd^xp{x,)f{x,)y Afjx,) n 

J d^x{dkf{xi)y d^xp{xi)f{xi) J d^x{dkf{xi)y^' 

This derivative vanishes for /(x*) = fmaxi^i) given by 

fmaxiXi) = CA~^p{Xi) = C'U{Xi), (40) 

where C and C' are inessential constants, which do not enter into the func¬ 
tional. It is easy to prove that the second variation of this functional is 
negative on the fmaxi^i), so this function provides the absolute maximum 
for the functional and due to the equations (34-36) its value coincides with 
the l.h.s. of (38). However, the function fmaxi^i) does not belong to the 
Hilbert space 11^2^ because the potential U{xi) has the asymptotic behaviour 
at infinity and therefore is not square integrable. 

The integral J which enters into inequality (38) could be written in the 
following form 

J = y d^xp{xi)fNixi), (41) 

where we denoted as fN^Xi) the normalized function f{xi): 


fN{Xi) 


f{Xi) 

f d^x{dkf{xi)y 


f{Xi) 

114/(3;*) II 2 ' 


(42) 


In the Hilbert space W 2 the Hoelder’s inequality holds true: 

I J d^xf{xi)g{xi)\ < ||/||p||^||p', ^ ^ ^ ^ 

as well, as the remarkable Ladyjenskaya’s [6] inequality: 

||/|U<(48)‘/«||4/I|2 (44) 
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Here we use the standard notations: 

II/lip = (/ d^xf{xify^\ 

\\9kf\\p = [j d^x\dkJ{xi)Y’y'^ 

From these two inequalities we obtain the following bound for the integral J: 

\j\ < (48)‘''»||p||«/5 (45) 

The facts mentioned above could possibly support the statement that the 
{—Estatic) is bounded from below by appropriate norm of the density. In¬ 
deed, the absolute maximum of the funcional \J\ should be bigger when its 
maximum in a restricted space like W 2 , given by (45). However we can not 
present the rigorous proof of this statement. 

One of the other general property of a static conhgurations of gas or flnid 
is the existence of the topological charge — ’’helicity” (or Hopf invariant), 
which explicit form is 

Q — J d^xv{xi)rotv{xi). (46) 

This object is not only the integral of motion of the equations (16) but it also 
is the central element of the algebra of Poisson brackets of {v{xi), p{xi)) [1]. 
The existence of such an object brings additional argument for the stability 
of the solitons. The role of ’’helicity” in the case of the solitons in plasma 
was pointed ont in [7],[8]. Moreover, in [9] it was shown that there exists 
a remarkable inequality which bounds the energy of plasma solitons from 
below by The derivation of snch inequality for our case (and in general 
for flnid solitons ) is highly desirable and we going to consider this question 
in the future publications. 

4 On the Possible Structure of the Static So¬ 
lutions 

The equations (32) which dehne the static conhgurations of gravitating gas 
are 3-dimensional nonlinear partial differential equations and the probability 
to hnd an analytic solution is very low. The only case where a class of soln- 
tions was fonnd in a similar sitnation is the t’Hooft-Polaykov monopole, bnt 
there the reqnirement of the spherical symmetry simplihed essentially the 
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problem. In our case we can not expect the spherically symmetric solution 
because the continuity equation requires the trajectories of the particles be 
closed in order to provide the static conhguration for v{xi),p{xi). The sim¬ 
plest and most symmetric conhguration we could expect for our case is the 
toroidal structure, where the density is concentrated in the vicinity of the 
axis of the toroid, while the held of velocity is tangential to the embedded 
one into the other toroidal surfaces. It is not the hrst time the toroidal-shape 
soliton appeared in the context of the theory of continuous media. Since the 
pioneer works of Lord Kelven in XIX century to the present time it was 
studied by many scientists both mathematician and physicists and recently 
the interest to the subject was again attracted by the works of Faddeev and 
Niemi [7],[8] 

It is clear that in the case of attractive gravitational interaction the par¬ 
ticles, which constitute the gas move around the region with bigger density 
(like planets around the Sun) so there should exists a collective motion in the 
approximation when the radius of torus tends to inhnity and we can speak 
about cylindrical rather when toroidal conhguration. Indeed, let us consider 
this axially symmetric tornado-like solution of (32). For that we shall write 
the Ansazt: 


v{xi) = v(r)(-p^,0), 

p{xi) = p(r), (47) 

where r = ^/x^ + y'^. The second and the third equations (32) are satished 
by (47), while the hrst gives 

^drv{rY=An—p{r). (48) 

r m 

This solution shows that the density stays an arbitrary function (what is ex¬ 
pected for the partial diherential equations) and the velocity grows with 
radius up to its asymptotic value. The later is the consequence of the 
approximation— here we actually have 2-dimensional potential logr instead 
of Coulomb K This example is to demonstrate that the particles which 
constitute the gravitating gas in their collective motion form the localized 
object. 

The traditional way to tackle 3-dimensional gas or huid is to introduce for 
the velocity Clebsh parametrization suggested in [10] and recently discussed 
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in[11], In general case this parametrization has the following form: 


v{xi) = f{xi)dg{xi) + dh{xi), (49) 

where f{xi),g{xi),h{xi) are scalar functions. For the toroidal solution with 
z - as the axis of symmetry we shall assume that the density p{xi) does 
not depend upon the azimuth angle (f) and Clebsh parametrization takes the 
form: 

v{xi) = cosa{r, z)d/3{r, z) + kd(j){xi), (50) 

where r, 0 are cylindrical coordinates, 


x2 _|_ y2 

Mxi) = arctang— 

X 

In fluid dynamics the parametrization (50) has 
interpretation which we shall discuss elsewhere, 
for this case has the following form 

q = f kd(f){xi) A dcosa{r, z) A d/3{r, z) 

= 27rk f dcosa(r, z) A d/3(r, z) (51) 

The function /3(r, z) in (50) should has a singularity at the point (tc, 0) in the 
semiplane (r, z) where Tc is the radius of axis line of the torus. When the point 
(r, z) goes around (vc, 0), the function /3(r, z) increases on 27rh The gradient 
of such function will be the vector tangential to the toroidal surfaces. The 
addition of the gradient of azimuthal angle </) makes the velocity (50) winding 
also in the azimuth direction. Such a behaviour of the held of velocity leads 
to the nontrivial Hopf invariant, which characterizes the homotopy classes 
7 r3(S'^). Described this way the Clebsh parameters, together with the density 
p(xi) could be found by numerical integration of the equations (32). 

The construction we presented and discussed above would be incomplete 
and too academic without an example of the galaxy which indeed may demon¬ 
strate the toroidal structure. Fortunately such galaxy does exists. It was 
discovered back in 1950 by Art Hoag and recently a very good picture was 
obtained by Hubble telescope (see hgure 1). On this picture it is clearly seen 

^Compare this with the similar parametrization of the dynamical variables in [8] for 
the case of plasma. 



an interesting mechanical 
The ’’helicity” functional 
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the toroidal structure and moreover the motion of the stars seem to be very 
close to what we should expect for the configuration with nontrivial ’’helic- 
ity” functional. In the gallery of the galaxies which is available on the sites 
http://www.astronomy.com or http://hubblesite.ogr we can hnd some other 
examples with the form more or less close to the torus, but the Hoag’s object 
is the best of all. these examples show that the model we considered might 
be not so trivial and could explain more sophisticated structures like spiral 
galaxy. 
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